Upstanding Rashba spin in honeycomb lattices: Electrically reversible surface spin polarization 
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The spin-split states subject to Rashba spin-orbit coupling in two-dimensional systems have long been ac- 
cepted as pointing inplane and perpendicular to the corresponding wave vectors. This is in general true for 
free-electron model, but exceptions do exist elsewhere. Within the tight-binding model, we unveil the unusual 
upstanding behavior of those Rashba spins around K and K' points in honeycomb lattices. Our calculation (i) 
explains the recent experiment of the Tl/Si(l 1 1)-(1 x 1) surface alloy [K. Sakamoto et ah, Phys. Rev. Lett. 
102, 096805 (2009)], where abrupt upstanding spin states near K are observed, and (ii) predicts an electrically 
reversible out-of-plane surface spin polarization. 

PACS numbers: 73.20.At,73.63.-b,71.70.Ej 



The honeycomb lattice is one of the three types of 
two-dimensional regular tessellation — triangular, square, and 
hexagonal tilings. In solid-state physics, the honeycomb lat- 
tice is described by two staggered triangular sublattices. The 
identity of the comprising atoms of the two sublattices deter- 
mines if the honeycomb is monoatomic or diatomic, and their 
interlayer distance determines whether the honeycomb is flat 
or bilayer. Thus graphene and boron nitride are monoatomic 
and diatomic fiat honeycombs, respectively, 1 Bi(lll) bilayer 
surface 2 is a monoatomic bilayer honeycomb, and Tl/Si(l 11)- 
(lxl) surface alloy 3 is a diatomic bilayer honeycomb. 

Among these honeycomb lattices, graphene has been under 
the most intensive investigation due to its unusual Dirac-like 
electronic excitations. 4 When deposited on a substrate, the 
structural inversion symmetry perpendicular to the graphene 
plane is broken, and the Rashba spin-orbit coupling 5 hence 
emerges. Recent experimental measurement of Rashba spin 
splitting in graphene on Ni(l 1 1) substrate with 6 and without 7 
an intercalated Au monolayer eventually drew theorists' at- 
tention to the electronic structure, in the presence of Rashba 
coupling, of graphene, 8 i.e., near K and K' in monoatomic 
flat honeycomb. In other honeycomb systems, Bi(l 1 1) bilayer 
surface is believed to contain strong Rashba coupling 9 and ex- 
hibit interesting spin-Hall patterns, 10 and a more recent exper- 
imental effort even shows an unusual Rashba spin behavior at 
K point in TI/Si(l 1 1)-(1 X 1) surface alloy. 3 

In this Rapid Communication we present a unified tight- 
binding description to understand the Rashba effect in hon- 
eycomb lattices. To focus on the major effect brought by the 
Rashba coupling, we consider electron hopping up to the near- 
est neighbors. Particular attention will be paid to the spin con- 
figuration near K and K' points, which shows an abrupt up- 
standing Rashba spin behavior, in good agreement with Ref. 
3. Moreover, we show that the upstanding spins along K 
and K' directions are opposite. Accordingly, we propose an 
electrically reversible out-of-plane surface polarization, which 
will be numerically shown by Landauer-Keldysh formalism. 

Consider a honeycomb lattice constructed by primitive 
translation vectors ti = a(l/2, V3/2, 0) and t2 = 
a(— 1/2, v3/2, 0), and basis vectors di = and d2 = 
a(0, 1 / \/3, d z ), where a is the lattice constant. For flat honey- 



combs we have d z = 0, and d~ ^ describes a bilayer case. 
We begin with the 4x4 tight-binding Hamiltonian 10 of the 
Slater and Koster type, 111 



I f 12 ^22 



(1) 



with off-diagonal element 



n 12 = 



U(l 
-it R (l - 



f 2F) 

F + V3G) 



-it R {l -F- V3G) \ 
U(1 + 2F) J 



(2) 

where tn is the Rashba hopping strength, U = l z V ppa + 
(l — lT) VppTr, l z being the direction cosine of near- 
est neighbors (l z = for flat and l z ^ for 
bilayer), is the two-center interaction integral involv- 
ing p z atomic orbitals, and the compact functions are 
given by F = exp(— i^/3k y a/2) cos (k x a/2) and G = 
exp(—iV3k y a/2) sin (k x a/2) . The diagonal elements of Eq. 
( 1) are Ma = E pi I with I the 2 x 2 identity matrix and i = 1, 2. 
For monoatomic honeycombs, we have E v \ = E. 
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for diatomic honeycombs, E p i ^ E P 2- In the following we 
consider E p i = E p 2\ straightforward generalization to the di- 
atomic case will be shown later. 

Adopting the same trick of Rashba, 8 Eq. (1) can be reduced 
to a 2 x 2 Hamiltonian, 
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for sublattice A, which depends explicitly on its eigen- 
value E. The Schrodinger equation of Hamiltonian (3) is 
ELa (E) \ipA) = E\iPa)- By solving the characteristic equa- 
tion det(ELA (E) — E) = the eigenvalues of Eq. (3) can be 
written as 



Eau = E p + (J,E V , jj,, v = ±1 
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where 

will play an important role in the following derivation. The 
eigenvectors can be written as either of 

m = (Vd 2 + \El - h n \ 2 y 1/2 ( ) (7a) 

(7b) 

which are independent of //. Both Eqs. (7a) and (7b) are 
valid for carrying out the spin expectation (S) = (h/2)(a), 
a = (<t x , <r y ,a z ) being the Pauli matrix vector, except at the 
symmetry points T, M, K and K' . With careful treatment at 
those points, the spin direction (S)^ = (h/2)(ip'^ / \a\ipj^) 
based on Eq. (7) subject to the four eigenstate branches are 
shown in Figs. l(a)-l(c); the band structure according to Eqs. 
(2)-(6) is shown in Fig. 1(d). Clearly one can see an abrupt 
upstanding component (S z ) near K and K' points. For 
the v = +1 branch, we depict the spin configuration based 
on Eq. (7) in Fig. 2, where each arrow is determined by 
((<7 X ) M+ , (a v ) IJ/ +) and the color shading is by (a z ) The 
bright (dark) region of (a z ) w 1 ((a z ) w -1) around K (K f ) 
can be clearly seen in the main panel of Fig. 2, reflecting its 
inherent C3 symmetry. 
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FIG. 2: (Color online) Spin configuration for the v = +1 branche. 
Each arrow is given by ({cr 1 }^, (<j v )i_i+) and the color shading is 
determined by (<t z ),j+. 



To provide deeper understanding of Fig. 1, as well as the 
zoomed-in plots of Fig. 2, below we give a series of discussion 
of £?„„ and around those symmetry points. 

The f point. Assuming k = 5 = <5(cos0, sin0) with 
5a <C 1, one can obtain 
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FIG. 1: (Online color) The spin components (a) {S x )nv, (b) {S y )^ v , 
and (c) {S z )fj.u of energy eigenstates of a Rashba spin-orbit-coupled 
monoatomic honeycomb lattice with its band structure shown in (d). 
The parameters are set tnj \U\ = 0.12. The Brillouin zone is shown 
in the inset in (d). 



and (<7) ;t „(f + ft) w — z/(sin cj), — cos </>, 0) which is identi- 
cal to the free-electron case. Equation (8) is useful for deter- 
mining parameters E p , U, and t r by matching with the free- 
electron dispersion E (k) = Eq + ti 2 k 2 /2m* ± ak. The band 
offset Eq, the Rashba parameter a, and the curvature h 2 /2m* 
(or the effective mass to*), as well as the lattice constant a are 
experimentally measurable. In the usual t 2 R <C |f/| 2 case we 
have.Ep_= £ + |3C/| , \U\ = 2h 2 /m*a 2 ,andtR = 2a/V3a. 

The M point. Let M = (2w/a,0), i.e., M 2 in Fig. 1. 
Assuming k = M + 5 = (2n/a + S x ,5 y ) one obtains an 
anisotropically free-electron-like dispersion, 

E V {M + 6) « yW44 - M (J2 _ 3^) a 2 

+ ^t R ^S 2 x + 98 2 a, (9) 

which indicates that M is a saddle point as one can see from 
the different sign of 5 2 and S 2 in the second term. This 
means that the effective mass of the electron at state near M 
has opposite sign when going along and perpendicular to the 
Brillouin-zone boundary. In addition, the latter has an effec- 
tive mass 3 times lighter and an effective Rashba parameter 
three times stronger than the former. Equation (9) therefore 
explains why we have band shape near f identical with that 
near M2 but not Mi [see Fig. 1(d)]. 

The anisotropy at M also reveals in the corresponding spin 
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where we keep terms up to first order in tfi/|J7|. Clearly from 
Eq. (10) the z component saturates to 2tn/\U\ along the Bril- 
louin boundary near M but vanishes when going perpendicu- 
lar to the boundary [see Fig. 1(c)]. 

The K and K' points. Let K' = (gi — g2)/3. Assuming 
k = K' + 5 = (47r/3a + S x , + 6 y ) we obtain 
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in agreement with Ref. 8. The spin direction near K', up to 
second order in (6a) , is given by 
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which shows at K' we have (a z ) M ± (i^'j = Tl- Around 
^ the dispersion is identical to Eq. (11), and the spin con- 
figuration has a reversed helicity, i.e., opposite out-of-plane 
component, 



(a*) ll „(K + 6) = -(a*) lxv (K' + 6), 



(13) 



but unchanged inplane component. 

So far the discussion is basically for sublattice A since 
we have obtained reduced Hamiltonian (3) by expressing the 
wave function of sublattice B, ips, in terms of that of sublat- 
tice A, ipA- We could have as well expressed ipA in terms 
of ipB', me resulting reduced Hamiltonian then would be 
Mb (E) = Epl + H f 12 Hi 2 / (E - E p ), leading to identical 
dispersion, identical inplane spin direction, but opposite (a z ) 
component, 



(14) 



In addition, we have so far focused on the monoatomic hon- 
eycomb lattice. For different sublattice atoms A and B, we 
have E p i ^ E p2 . The energy dispersion can be straightfor- 



wardly derived as e pv (k) = E p + \ie v with e v 



E 2 



E p = (E p i+E p2 )/2, A p = (E p i-E p2 )/2, and E v given by 
Eq. (5). Noticing from Eqs. (8), (9), and (11) that around the 
symmetry points E+ ^ E- occurs only at K and K' , we ex- 
pect the modification to the band structure due to A p appears 
the most salient at K and K' . Indeed, with E v (\t = K or 
K') = M R (1 + i/)/2a gap e+_ - = 2e_ = 2 |A P | is 
opened. At other symmetry points, e p ± remain degenerate. 
Despite the opened gap, however, the eigenvectors surpris- 
ingly remain the same as Eq. (7), and the upstanding Rashba 
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FIG. 3: (Color online) Energy dispersion and {a z ) l±v (indicated by 
arrows) for (a) monoatomic and (b) diatomic honeycomb lattices. In 
(a) the black dashed lin es are the approximated dispersions given 
by Eq. (11); E M = + 4t 2 R . In (b) we set A p = 0.1 \U\; 



em = + U 2 + 4t 2 R and e K = \J A 2 + (3t R ) 2 . The inset 

shows a six-terminal bilayer honeycomb channel with black (gray) 
dots the upper (lower) sublattice. 



spin around K and K' described previously is therefore un- 
changed. 

To show the opened gap with unchanged spin direction, 
we plot the energy dispersion curves with (a z ) fllJ in Fig. 3 
for monoatomic honeycomb near K' [panel (a)] and diatomic 
honeycomb near K [panel (b)]. Figure 3(b) readily explains 
the recent spin-resolved angle-resolved photoelectron spec- 
troscopy data for the Tl/Si( 1 1 1 )-( 1 x 1) surface alloy. 3 In such 
surface alloy, 12,13 the Tl coverage on the Si(lll) substrate is 
one monolayer, and the topmost Tl (sublattice A) and Si (sub- 
lattice B) layers form a diatomic honeycomb lattice of the 
bilayer type (d z < 0). The \l — +1 branch in Fig. 3(b) thus 
resembles the band feature near K reported in Ref. 3. It is 
important to note, however, that for flat honeycombs such as 
graphene, both sublattices A and B will be silmultaneously 
measured, leading to vanishing (S z ), in view of Eq. (14). 

We have shown that the upstanding Rashba spin around 
K and K' points is a fundamental property of the spin 
configuration in honeycomb lattices, whether flat or bilayer, 
monoatomic or diatomic. Next we illustrate how striking this 
property can be. For simplicity, let us consider a six-terminal 
channel made of monoatomic bilayer honeycomb. See the in- 
set in Fig. 3. Assume that the transport of this six -terminal 
device is supported by the surface states so that electrons are 
only allowed to hop within the surface bilayer, i.e., the hon- 
eycomb lattice. Let the Fermi energy Ep lie just at E p . The 
situation is like Fig. 3(a). Recalling Eq. (14) we expect that 
when driving the electrons along K', the transport states with 
Ep = Ep are +z spin polarized on sublattice A but — z spin 
polarized on sublattice B. Recalling further Eq. (13), the sur- 
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FIG. 4: (Online color) Surface spin polarization in a 6-terminal bi- 
layer honeycomb channel with (a) A"'biasing, (b) and (c) K bias- 
ing, and (d) M biasing. Red/gray (green/light gray) dots represent a 
(<j z ) > 0((a z ) < 0) local spin density with the dot size proportional 
to |(cr z )|. (ct 2 ) and (erll) as functions of the bias eVo are shown in 
(e) for KIK 1 biasing and in (f) for M biasing. 



face spin polarization (assumed to be contributed by sublattice 
A only) is out-of-plane and can be electrically controlled by 
either K or K' biasing. 

To visualize this idea, let us calculate for the six-terminal 
bilayer honeycomb channel the local spin densities by em- 
ploying the Landauer-Keldysh 1415 formalism, subject to 
Hamiltonian, U = E p J2 n c l c « + Y,(nm) c \n[ ul + x 
d rn „) z ]c„, with cjj (c„) the creation (annihilation) operator of 
the electron on site n and d mn the unit vector pointing from 
site n to m. The total number of lattice sites in the honey- 
comb channel is 480. For clarity we will plot (a z ) only on 



the surface (sublattice A). A positive (negative) (a z ) on each 
site will be denoted by a red (green) dot, with the dot size 
proportional to the magnitude of (<r z ). The applied potential 
energy of ±eVb/2 will be denoted as "±," and eVb = as 
"0" on each lead. The parameters E p = 0, U = — 1 eV, and 
£r/ \U\ = 0.12 are identical to those used in previous figures, 
and are within a realistic range. 16 The bias is eVo = 2 mcV. 

First we drive the electrons from left to right, corresponding 
to K' . As expected, we have a positive average of surface spin 
polarization (a z ) > 0, as shown in Fig. 4(a). When rotating 
the bias direction by 60° counterclockwise, the surface polar- 
ization becomes (a z ) < as a consequence of Eq. (13) [see 
Fig. 4(b)]. Reversing the bias of Fig. 4(a) also switches K' 
states to K, leading to (a z ) < 0, as shown in Fig. 4(c). In Fig. 
4(d) we drive the electrons from bottom to top, corresponding 
to the M direction. The spin density distribution becomes 
completely different and satisfies the intrinsic spin-Hall sym- 
metry, which yields (<j z ) = 0. Figures 4(a) and 4(c) imply that 
the out-of-plane surface spin polarization can be flipped sim- 
ply by reversing the bias, which is also a direct consequence 
of time-reversal operation. To show this electrical control of 
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surface spin, we plot (a z ) and (cry) = {{a x ) + {a y ) ) 1 / 2 
as a function of the bias eVb in Figs. 4(e) and 4(f) for K' IK 
and M biasing, respectively. Within the low bias regime [see 
the inset between Figs. 4(e) and 4(f)], (a z ) for K' IK biasing 
grows with eVb linearly, while other components are either 
vanishing or relatively small. 

In conclusion, we have presented a unified tight-binding de- 
scription to understand the Rashba effect in graphene, as well 
as bilayer surfaces and surface alloys of the honeycomb struc- 
ture. Our results explain the recently observed abrupt upstand- 
ing Rashba spin in Tl/Si(l 1 1)-(1 x 1) surface alloy around K, 
and predict an electrically reversible out-of-plane surface spin 
polarization, which may serve as a storage mechanism for fu- 
ture spintronic devices. 
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